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ABSTRACT 

Context. Realistic stellar models are essential to the forward modelling approach in asteroseismology. For practicality however, certain 
model assumptions are also required. For example, in the case of subdwarf B stars, one usually starts with zero-age horizontal branch 
structures without following the progenitor evolution. 

Aims. We analyse the effects of common assumptions in subdwarf B models on the g-mode pulsational properties. We investigate 
if and how the pulsation periods are affected by the H-profile in the core-envelope transition zone. Furthermore, the effects of C- 
production and convective mixing during the core helium flash are evaluated. Finally, we reanalyse the effects of stellar opacities on 
the mode excitation in subdwarf B stars. 

Methods. We computed detailed stellar evolutionary models of subdwarf B stars, and their non-adiabatic pulsational properties. 
Atomic diffusion of H and He is included consistently during the evolution calculations. The number fractions of Fe and Ni are 
gradually increased by up to a factor of 10 around log T = 5.3. This is necessary for mode excitation and to approximate the resulting 
effects of radiative levitation. We performed a pulsational stability analysis on a grid of subdwarf B models constructed with OPAL 
and OP opacities. 

Results. We find that helium settling causes a shift in the theoretical blue edge of the g-mode instability domain to higher effective 
temperatures. This results in a closer match to the observed instability strip of long-period sdB pulsators, particularly for / < 3 modes. 
We show further that the g-mode spectrum is extremely sensitive to the H-profile in the core-envelope transition zone. If atomic 
diffusion is efficient, details of the initial shape of the profile become less important in the course of evolution. Diffusion broadens the 
chemical gradients, and results in less effective mode trapping and different pulsation periods. Furthermore, we report on the possible 
consequences of the He-flash for the g-modes. The outer edge of a flash-induced convective region introduces an additional chemical 
transition in the stellar models, and the corresponding spike in the Briint-Vaisala frequency produces a complicated mode trapping 
signature in the period spacings. 

Key words, subdwarfs - stars: evolution - stars: oscillation - methods: numerical 



1. Introduction 

Hot B-type subdwarfs are identified as core He-burning stars sur- 
rounded by a very thin H-envelope (Heber 1986). This places 
them at the blue extension of the horizontal branch in the 
Hertzsprung-Russell diagram. Hence, they are also referred to as 
extreme horizontal branch (EHB) stars. Subdwarf B (sdB) stars 
are ubiquitous in our Galaxy, where they dominate the popula- 
tion of faint blue objects at high galactic latitudes ( |Green et al. 
|1986 1. They are also believed to be responsible for the ultravi- 
olet excess or 'UV upturn' in giant elliptical galaxies ( |Brown| 
|et al.|199T |Yi et al.|1997| l. A thorough review of hot subdwarfs 
is given by Heber (2009 1. 

The future evolution of sdB stars is straightforward and 
undisputed. After He in the core is exhausted, a short phase of 
He-shell burning follows, and the star may be identified as a hot- 
ter subdwarf of O-type (sdO). The H-envelope is too thin to sus- 
tain H-shell burning, so the star will end as a C-O white dwarf 
without evolving through the asymptotic giant branch phase. 
The formation of sdBs on the other hand, remains a much de- 
bated topic. Many formation channels have been proposed such 
as enhanced mass loss on the RGB ( |D'Cruz et al.|19 96), mass 
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loss through binary interaction (Mengel et al. 1976| l, and the 
mergers of two He white dwarfs (Webbink 1984). The relative 
importance of different formation scenarios has been evaluated 
by binary population synthesis studies ( Han et al.|2002 2003 i. 
While these studies are valuable, it should not be forgotten that 
they are dependent on parametric descriptions of binary inter- 
action. Detailed mass determinations of subdwarf B stars could 
provide important constraints on the binary interaction mecha- 
nism. Accurate astrophysical mass determinations are only pos- 
sible in special circumstances such as in eclipsing binary sys- 
tems. Asteroseismology provides an alternative method, since 
it allows a detailed study of the interior of non-radially pulsat- 
ing stars. The consistency of both approaches has been achieved 
for PG 1336-018 (see |Vuckovic et al.|2007| |Hu et al |2007l and 
Charpinet et al.|2008[ ). 

Subdwarf B stars exhibit a variety of pulsations. The first 
sdB pulsator, EC 14026-2647, was observed by Kilk enny et al.| 
( 1997 ) to pulsate in multiple short-period modes. This prototype 
represents the variable class V361 Hya stars which now includes 
42 rapid sdB pulsators, with periods in the range 80-600 s. 
The short -period modes have been interpreted as low radial or- 
der, low spherical degree /?-modes (Charpi net et al.|19 97). The 
V361 Hya stars are amongst the hottest sdB stars with effective 
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temperatures between 28 000 K and 35 000 K and surface grav- 
ities 5.2 < logg < 6.1. At the cooler end of the EHB, 31 sdB 
pulsators with periods between 30 min and two hours have been 
discovered by Gre en et aT[ ([2003 ). This variable class is named 
V1093 Her but is also often referred to as PG 1716+426 stars af- 
ter the class prototype. The long periods suggest that these stars 
pulsate in high radial order g-modes. Especially interesting are 
the sdB pulsators that exhibit both p- and g-mode pulsations. 
Three of these so-called hybrid pulsators have been found at the 
intersection of the V361 Hya and VI 093 Her stars (Oreiro et al. 



[20051 |Baran et al.|2005[|Schuh et al.|2006l|Lutz et al.| 2008 ) 

The opacity mechanism operating in the Fe opacity bump 
around logT = 5.3 has been successfully used to explain the 
excitation of p- as well as the g-modes (Charpinet et al. 1996 



Font aine et al.|2003 i. This opacity bump is caused by Fe accu 



mulation owing to the competing diffusion processes of gravita- 
tional settling and radiative levitation. While seismic mass de- 
terminations have been achieved for a dozen V361 Hya stars us- 
ing static envelope models (see Fonta ine et al.| [2008 and refer- 
ences therein), this is not yet the case for V1093Her stars. The 
reason is twofold. First of all, the g-modes have lower ampli- 
tudes and longer pulsation periods and thus require higher pre- 
cision observations and longer observing runs to detect the fre- 
quencies. Secondly, envelope models may suffice to describe the 
V361 Hya stars' shallow p-modes that probe only the outer lay- 
ers. However, they are not suitable for modelling g-modes that 
penetrate to deeper stellar regions. Stellar models with detailed 
information about the He-rich core are required to predict reli- 
able g-mode periods. 

Another problem for the VI 093 Her stars is the discrep- 
ancy between the observed and theoretically predicted instabil- 
ity strip. The first theoretical models that show unstable g-mode 
pulsations ( Fontaine et al.|2003 1 have much lower effective tem- 
peratures than observed with a discrepancy of ~ 5000 K. Jeffery 
& Saio (2006b ) managed to place the blue edge of the instability 



strip within ~1000 K of the observed blue edge. They did this 
by using OP opacities (|Badnell et al.||2005] l rather than OPAL 
opacities (Iglesias & Rogers 19961, and by enhancing the Ni 
abundance in the envelope in addition to Fe. Thus, their work 
suggests that sophisticated models including up-to-date opaci- 
ties and diffusion of other Fe-group elements are required to ad- 
dress the instability issue properly. 

An exploratory study of how the pulsational properties of 
sdB stars are related to their internal structure has been carried 
out by Charpinet et al. (2000). They showed that the He-H tran- 
sition zone between the He-rich core and H-rich envelope is re- 
sponsible for trapping and confining g-modes in sdB stars. This 
is similar to the well-known mode trapping phenomenon in com- 
positionally stratified white dwarfs (Wing eFet al.||198l} . Mode 
trapping results in deviations from the asymptotic constant pe- 
riod spacing, making the g-modes sensitive to the mass of the 
H-envelope. Furthermore, the C-O/He-transition region between 
the convective and the radiative part of the core will also influ- 
ence the g-mode period distribution. This may be a new interest- 
ing way to follow the He-core evolution, which would also have 
implications in a broader context for horizontal branch stars. 

Here, we shall explore the sensitivity of the g-modes to cer- 
tain assumptions about the He-flash. This energetic event charac- 
terizes the start of He-ignition in degenerate cores of low-mass 
stars (<2 M ). Most stellar evolution codes encounter numeri- 
cal difficulties in calculating this phase, and even in successful 
cases uncertainties remain because of the approximative treat- 
ment of convection. Three-dimensional hydrodynamic calcula- 
tions are necessary to investigate this enigmatic phase of stel- 



lar evolutio n (e.g., |Deupree|19 96; Dearborn et al.|200 6; Mocak 
et al.||200"8 ), but these are expensive in computing time. When 
many models are needed, as in forward modelling, it is more 
practical to construct post-flash models without following the 
previous evolution in detail. One should then use the outcome of 
detailed studies to ensure that the stellar models are as realistic 
as possible. Most simulations indicate that the initial flash starts 
off-centre, followed by several mini-flashes moving towards the 
centre. While it is generally found that most of the inner region 
will be convectively mixed, there are uncertainties in the outer 
extent of the convective region and the amount of He-burning 
during the He-flash varying from 3% to 7% in mass fraction 
( |Piersanti et al.|2004||S"erenelli & Weiss|2005| ). The edge of the 
He-flash-induced convective region will leave a chemical com- 
position gradient, which we believe can produce additional mode 
trapping features. 

We start with constructing zero age extreme horizontal 
branch (ZAEHB) models in Sect. [2] We make an analytical fit to 
the H-profile in the core-envelope transition region and examine 
its influence on the g-mode periods (Sect. |2~2] i. The effects of un- 
certainties in the He-flash treatment, i.e. convective mixing and 
C-production, are explored in Sect. 2.3 In Sect. [3] we perform a 
stability analysis on a grid of evolutionary sdB structures during 
the core He-burning phase. We show the impact of the opacities 
and H/He-diffusion on the instability strip. The conclusions are 
summarized in Sect. [4] 

2. Zero age extreme horizontal branch models 

2.1. Approach 

The stellar models were calculated using the stellar evolution 
code STARS ( |Eggleton|197T]|Han et al.|T994}|Pols et al.|1995| ). 
The details of the input physics used in this work are described 
in Hu et al. ( 2009 ). By default, we constructed the stellar models 



with OPAL opacities (except in Sect. 3.3 i. The pulsational prop 



erties were calculated with the adiabatic oscillation code OSC 
( |Scuflaire et al. 2008|) in combination with the non-adiabatic 



Dupret (2001) 



code MAD by 

Our procedure of constructing ZAEHB models is as follows. 
First, we take the stellar structure model of a higher mass (2.25 
M Q ) model that just ignited helium non-degenerately. We then 
impose an artificial mass loss (10~ 6 M* yr~') until the remaining 
mass is of our choice. Finally, we allow the star to evolve, while 
slowly changing the chemical abundances to the values we wish 
(see Sect. |2.2| i, and after switching off composition changes pro- 
duced by nuclear reactions and convective mixing. We ensure 
that the new ZAEHB model is fully converged by checking the 
increase in the timestep that the evolution code gives for con- 
verged models. 

2.2. The He-H chemical transition profile 

Commonly, one constructs ZA(E)HB models without follow- 
ing the previous evolution. However, special care must be taken 
when reconstructing the stellar chemical abundances. In partic- 
ular, the shape of the H-profile in the core-envelope transition 
zone is very important to g-modes because of the associated 
mode trapping. Therefore, we first evolved a 1.00 M ZAMS 
model to the tip of the red giant branch and made a careful fit to 
the H-profile of the red giant (see top panel of Fig. [Tj(. We sug- 
gest a formula that fits the H-profile very well using two third 
degree polynomials: 



X(M r ) = 



for M r < Mi, 



(1) 
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(A(M r -Mi)Y 
(B(M r - M 2 )f + X surf 

^surf 



for Mi < M r < Mii nk , 

for Miink < M r < M 2 , 
for M r > M 2 . 



We found that the values M\ = 0.46575 M , M 2 = 0.4661 M G , 
M Iink = Mi + 1.5 x 10~ 4 M Q , X surf = 0.688, and A = 4400 
provide the best-fit profile. The coefficient B is then determined 
by assuming that the curve is continuous at M r = Mi + Mn^. 

We compared our approximation to a sinusoidal profile as 
proposed by Sweig art & Gross] ( 1976) and used in the ZAHB 
models of Jeff ery & Saio| ( |2006a| l. The H abundance in the tran- 
sition zone according their approximation is given by: 



X(M r ) = 0.5X surf 



1 - cos 



\ M 2 -Mj 



for Mi < M r < M 2 .(2) 



In the top panel of Fig. [T] we show both approximations, and 
find that our fit provides a more accurate representation of the 
red giant's H-profile. 

We constructed three ZAEHB models of total mass 0.470 
M Q , i.e., with a H-rich envelope of 9 x 10~ 4 M above the He-H 
transition zone. In model (i), the abundances were taken from 
the red giant model throughout the star. In models (ii) and (iii), 
the H-profile is given by our suggested Eq. ([l} and by the sinu- 
soidal Eq. (pi, respectively. In the latter two cases, we simply 
assume that Z = 0.02 and there is a |Grevesse & Noels] (1993) 
metal mixture throughout the star. In model (i), there are minor 
effects on the metallicity because of CNO cycling, atomic diffu- 
sion, and convective mixing in the previous evolution. We note 
that we have not yet accounted for flash-induced C-production 
or mixing. This is to ensure that our comparison is meaningful 
here. The effects of the He-flash will be evaluated separately in 
Sect. |3] 

It is well known that in the asymptotic limit the g-mode pe- 
riods are equally spaced in the radial order n (Tassoul 1980). 
However, sharp features in the Briint-Vaisala frequency cause 
mode trapping, which manifests itself in deviations from the uni- 
form period spacing. For an in-depth treatment of mode trapping 
in white dwarfs, sdB stars, and main-sequence stars, we refer the 
interes ted r eader to e.g., Brassard et al.| ( |1992| i, |Charpinet et al.| 
(2000), and Miglio et al. (2008}, respectively. Here, it is impor- 
tant to realize that the period spacing between two consecutive 
modes in radial order (APi = P; i/7+ i - P;,„) can be described by 
an oscillatory component superposed on the asymptotic constant 
period spacing. The periodicity An of the oscillatory component 
is related to the location of the chemical gradient, while its am- 
plitude is determined by the steepness of the gradient. Trapped 
modes are then located at the minima in the AP versus P dia- 
gram. From the asymptotic theory (Tassoul 1980), one can de- 
rive the approximative relations 



An = 



Fie 

n n 



l,i+&n 



n H 



V/(/ + i)' 



where 



n H = 2k 



and n = 2k 2 



(3) 



(4) 



The Briint-Vaisala (BV) frequency N is given by 
2 _ GM r j 1 dh\P dlnp\ 



N = 



dr 



dr 



The integration boundary ru indicates the location of the He-H 
transition zone, r c the edge of the convective core (= for a 



red giant H-profile 
our fit 
sinusoidal fit 
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Fig. 1. The effect of the H-profile on mode trapping. Upper 
panel: the H-profile in the core-envelope transition zone. We 
compare the H-profile of a red giant model (thick solid line) 
to our polynomial Eq. ([T]( (blue dashed line) and to the sinu- 
soidal Eq. Q (thin solid line). The three lower panels give the 
period spacing between modes of increasing radial order n as a 
function of the pulsation period for the modes of spherical de- 
gree / = 1, 2, 3. We indicated with the horizontal dotted lines the 
asymptotic constant period spacing. 



ZAEHB model), and R is the total stellar radius. The asymptotic 
constant period spacing is given by YIq/ y/l(l + 1). 
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Figure [T] shows the sensitivity of the g-modes to the exact 
shape of the H-profile. We plotted the period spacing AP against 
the pulsation period P, for the modes I = 1,2,3. We see that 
our polynomial fit to the H-profile reproduces the g-mode spec- 
trum of model (i) very well. In contrast, the sinusoidal H-profile 
leads to a significantly different mode trapping signature. The 
deviations become more distinct at higher spherical degree I and 
higher radial order n. Another interesting feature can be seen 
especially for I — 3 in the AP versus P diagrams, where there 
are two oscillatory components in the period spacing. We see a 
short periodicity in the radial order of An\ = 3, modulated by a 
long periodicity of A«2 = 28. This is at first instance unexpected, 
since we only have one chemical transition zone in our ZAEHB 
models. 

Numerical integration of Eq. Q for our model (ii) gives 
An = 3, which accounts for the short periodicity. The long pe- 
riodicity, however, cannot be explained by the above relations. 
To understand this effect, one should realize that trapping oc- 
curs for a mode if the nodes of its eigenfunctions are suitably lo- 
cated with respect to the BV-discontinuity. In other words, only 
modes with a certain phase with respect to the discontinuity are 
trapped, which gives rise to the periodicity of Eq. ([3]). In reality, 
the BV-discontinuity of course has a finite width, which was not 
accounted for in the derivation of Eq. Q. Therefore, depend- 
ing on the change of phase within the BV-discontinuity, certain 
modes can be trapped more efficiently than others, leading to the 
second periodicity. Again from the asymptotic theory of Tassoul 
( 1980), but now taking into account the finite width of the BV- 
spike, it can be shown that approximately (Hu et al., in prepara- 
tion), 

A» 2 = ^ with Il m = 27C 1 (£^drj\ (5) 

where r\ is the radius at the inner boundary of the BV-spike and 
r2 is at the outer boundary. Numerical integration of Eq. <[5j for 
our model (ii) indeed results in A«2 = 28, showing that the ap- 
proximate relations hold up very well. 

To clarify this, it is helpful to examine the displacement 
eigenfunctions, f r (radial) and (horizontal). We find that a 
mode is trapped when it has a node in £/, at the bottom of the 
He-H transition zone directly followed by a node in £ r . This is 
consistent with the findings of Brassa rd" et al.| ([T992) for white 
dwarfs. In Fig. [2] we plot for model (ii) the displacements for 
three trapped modes of radial order n = 36, n = 45, and 
n — 36 + A«2 = 64. These have a small, large, and again small 
amplitude in AP, respectively (see Fig. [TJ. We see that mode 
n = 36 also has a node in £/, at the top of the transition zone, 
while mode n = 45 has a node there in £ r . We find this to be gen- 
erally true and conclude that modes are more efficiently trapped 
when the first node encountered in the He-H transition zone cor- 
responds to £a and the last node to tj r . After A«2 modes, the nodes 
are again suitably placed. This can be seen by comparing n = 36 
and n = 64: the location of the nodes of with respect to the 
boundaries is approximately the same, while for n - 64, there is 
one additional node in-between. The interpretation of this inter- 
esting phenomenon will the subject of a forthcoming paper (Hu 
et al., in preparation). 




-0.2 1 1 1 1 1 

0.2 0.4 0.6 0.8 

r/R. 



Fig. 2. Normalized displacement eigenfunctions £ r (dashed line) 
and £/, (solid line) as a function of the normalized stellar radius 
for ZAEHB model (ii). The upper panel is for mode (I, n) = 
(3, 36), the middle panel for (/, n) = (3, 45), and the lower panel 
corresponds to (I, n) = (3,64) (see Fig. [T] bottom panel). We 
zoom into the He-H transition zone to show the behaviour of the 
eigenfunctions in the g-mode cavity. In the close up, we also plot 
the H-profile (dotted line) to indicate the extent of the transition 
zone. 



2.3. The helium core flash 

In previous models we did not account for the inner core abun- 
dances being in reality affected by nuclear burning and convec- 
tive mixing during the He-flash. From now on, we treat the He- 



flash in a 'standard' way by converting 5% of the He-mass frac- 
tion into C and assume that the He-flash-induced convection re- 
gion extends to 90% of the He-rich core (i.e., M r = 0.42 M Q ) 



as suggested by 3D hydrodynamical simulations of Dearborn 
|et al.| ( [2006] >. To test the influence of nuclear burning and convec- 
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Fig. 3. Effects of different assumptions about the He-flash on the g-mode spectrum of ZAEHB models. Left: the BV-frequency 
throughout the star in terms of the temperature. The dashed curves are shifted upwards with log 500 for visibility. Right: the period 
spacing between consecutive g-modes as a function of the pulsation period. In the upper panels the case of C-production during 
the He-flash of 5% throughout the entire He-core (solid line) is compared to the case of no C-production (dashed line). The middle 
panels show the results for ZAEHB models where convective mixing was assumed up to M r = 0.42 M (solid line) and M r = 0.45 
M Q (dashed line). In the lower panels the amount of C-production is varied from 3% (dashed line) to 7% (solid line). 



tive mixing during the He-flash, we also constructed post-flash 
ZAEHB models where i) C is increased by 3% and 7%, respec- 
tively, and ii) the core is convectively mixed with up to 97% 
(0.45 Mq) and 100% (0.46575 M ) of the core, respectively. 

Figure|3]shows the effects of the different assumptions about 
the He-flash on the g-mode spectrum. The upper panels indicate 
that increasing C in the entire He-core does not have a significant 
effect. There is only a minor shift to higher pulsation periods be- 
cause of a slight increase in the stellar radius, but the g-mode 
spectrum remains qualitatively the same. However, if the He- 
flash-induced C-production is not mixed throughout the entire 
core, there is an additional chemical gradient, and thus an addi- 
tional spike in the BV-frequency. This results in a complicated 
mode trapping pattern in the period spacings, that depends on the 
location of this spike, as shown in the middle panels of Fig. [3] In 



the lower panels, one can see that the strength of the BV-spike, 
and thus the effect of mode trapping, is sensitive to the amount 
of C produced during the He-flash. 



We emphasize that the ZAEHB models do not yet have a C- 
O convective core, and thus the He-H transition and the edge of 
the He-flash-induced convective region are mainly responsible 
for mode trapping in these models. The thin He-convection zone 
near the surface, which is visible in Fig.|3]as negative N 2 around 
log T = 4.7, has only a very small effect on the g-modes. 
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Fig. 4. Comparison of sdB models with (dashed line) and without (solid line) atomic diffusion at an EHB age of 5 x 10 7 years. The 
left panel shows the BV-frequency as a function of the stellar interior temperature. The dashed curve is shifted upwards with log 500 
to avoid overlapping. The three BV-spikes from left to right correspond to 1) the boundary of the convective core, 2) the edge of 
the He-flash-induced mixing region (small spike), and 3) the He-H transition zone. The middle panel shows the effect of diffusion 
on the H-profile. The right panel shows the period spacing versus the period, illustrating the mode trapping caused by the chemical 
transitions. 



3. Evolutionary extreme horizontal branch models 

3.1. Approach 

During the EHB evolution, Fe accumulates in the driving region 
(around log T = 5.3) due to radiative levitation (Charpinet et al. 
1997 Fontaine et al.|200"6 l. It is reasonable to expect that other 
Fe-group elements, such as Ni, are affected similarly by radia- 
tive levitation. Considering the work of |Jeffery~& Saio (2006b), 
we included both Fe and Ni enhancement in our sdB models. 
We built tables from OPAL opacities (Iglesias & Rogers 1996), 
where the number fractions of Fe and Ni are enhanced by a fac- 
tor / = 1, 2, 5, 10 relative to the |Grevesse & Noels| ( |1993| l metal 
mixture. Of course, the number fractions of the other metals are 
accordingly reduced. The OPAL radiative opacities are supple- 
mented by low -temperature opacities (Ferguson et al. 2005) and 
conductive opacities ( Cassisi et al.|2007 l. 

We increased the enhancement factor / during the EHB evo- 
lution from 1 to 10 with a Gaussian centred around log T = 5.3. 
Already after 10 7 yr, we have /(log T = 5.3) ~ 8 (see Fig. 1 
of Hu et al. 2008). The opacity is then calculated by interpo- 
lating between the opacity tables for different /. We note that 
Fe/Ni-accumulation is only accounted for in the opacities. This 
is reasonable since the mass fractions of Fe and Ni are low even 
after enhancement, and these elements are not involved in nu- 
clear reactions. Thus, the only relevant effect of the enhance- 
ment is in the stellar opacities. We follow the approach by |Hu| 
et al.| ( [2008] l, but with two noteworthy differences taking into ac- 
count the work of |Jeffery & Sa io (2006b). First of all, as we have 
already mentioned, Ni is enhanced in addition to Fe. Secondly, 
the mass fractions of the non-enhanced metals remain constant 
resulting in a slight increase in the total metallicity, whereas in 
our previous work, the other metals were reduced to keep the 
total metallicity constant. 



3.2. Atomic diffusion 



Our models include H/He-diffusion as described in |Hu et al7| 
(2009). In the presence of diffusion, steep composition gradi- 
ents are broadened leading to less efficient mode trapping. This 
can be seen in Fig. [4j where we show the same ZAEHB model 
evolved with and without atomic diffusion. Both models have an 
EHB age of 5 X 10 7 years. As expected, the amplitude of the 
oscillatory component in the period spacing is much smaller in 



the case of diffusion. Because of the lower amplitude, the ef- 
fect of A«2 (~ 30 in this case) is less significant, although still 
visible. The mode trapping pattern is at this stage even more 
complicated, because the sdB star has developed a convective 
core, which leads to another BV-discontinuity (see left panel of 
Fig.g). 



In Sect. 2.2 we detected the extreme sensitivity of the g- 
mode spectrum to the shape of the H-profile. This shows that 
special care must be taken when modelling the He-H transition. 
However, as atomic diffusion tends to wash away the differences 
in the initial chemical profiles, the significance of modelling the 
H-profile very precisely becomes less for evolved models. 

Another interesting effect of He-settling is the disappearance 
of the He-convection zone, as can be seen in the left panel of 
Fig. |4] The negative A^ 2 around log T = 5.3 corresponds to the 
Fe surface convection zone. We note that the Fe convection zone 
was not present in our ZAEHB models, since our models accu- 
mulate Fe during the evolution. It is not expected that the surface 
convection zones have a significant effect on the g-mode periods 
because of the small contributions of the surface layers to the 
weight function (see Fig. 9 of Charpin et et al.|20 00). 

3.3. Stability analysis of a grid of sdB models 

The previous models all have a mass of 0.470 M G . For a sta- 
bility analysis, we constructed models of varying mass so that 
the observed VI 093 Her logg - log T e s region is mapped. We 
computed sdB models with a fixed core mass defined as before 
at 0.46575 M Q , and envelope masses: 2.5 x 10~ 4 , 5.5 x 10~ 4 , 
1.25 x 10~ 3 , 2.25 x 10~ 3 , 4.25 x 10~ 3 , 6.25 x 10~ 3 , and 8.25 x 10~ 3 
M . To make the models as realistic as possible, we took the 
abundances from the red giant model at the tip of the RGB. The 
inner 0.42 M Q of the core is assumed to be mixed during the 
He-flash and C production is taken to be 5%. 

It is well-known that mode excitation is very sensitive to the 
stellar opacities. Thus, we also constructed opacity tables with 
OP opacities using OP server ( [Badnell et al.|2005[|Seaton|2005| l. 
For a fair comparison, we obtained OP tables for the same metal 
mixture as OPAL (i.e., |Grevesse & N oels 1993). Because the 
maximum temperature for OP opacities is logT = 8, we use 
OPAL opacities for He-burning regions that are enriched in C 
and O. 

Figure [5] shows the evolutionary tracks computed with OP 
opacities in the log g - log T s g diagram. The tracks for OPAL 
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overlap with these, thus the differences between OP and OPAL 
opacities have a negligible effect on the global evolution. The 
structure models for which we calculated non-adiabatic pulsa- 
tional properties, are represented by circles. The time interval 
between the pulsation models is 10 7 years. 



Table 1. Blue edge in terms of log r e jf of the g-mode instability 
strip for our sdB models. 





4.6 - 




4.8 ; 




5 - 


CM 








E 


5.2 - 


C3, 




CT 




O) 

_o 


5.4 - 




5.6 - 




5.8 - 



6 

4.55 




4.45 4.4 

log T eff (K) 



4.25 



Fig. 5. Evolutionary tracks of our sdB models in the logg - 
log T e g diagram. The tracks from high to low T e g correspond 
to sdB masses increasing from 0.466 to 0.474 M . It is assumed 
that 5% of C is produced during the He-flash and the inner region 
is convectively mixed up to M r = 0.42 M Q . The circles indicate 
the models used in the stability analysis. The dotted box shows 
the region in which VI 093 Her stars have been observed. 



Figure [6] shows the instability strip of our models obtained 
using OPAL opacities (upper panels) and OP opacities (mid- 
dle panels). In agreement with previous studies (Jeff ery & Saio| 
2006b] |Miglio et al.|2007) , we find that with OP opacities, un- 
stable modes are found in hotter stars than with OPAL opacities. 
This is related to differences in the strength and location of the 
Fe-group opacity bump as a function of the temperature. For OP 
opacities, this bump is slightly larger and occurs at higher tem- 
peratures (see Seaton & Badnell 2004). The thermal relaxation 
time (Tth = CyTIS.mlL,) in the driving region, i.e., around the 
opacity bump, is then large enough for efficient driving of the 
g-modes. This is because, first of all, the opacity driving mecha- 
nism works most efficiently on modes with periods comparable 
to T t h. Secondly, only high-order g-modes (with long periods) 
can be excited, for low-order g-modes are affected by significant 
radiative damping. The same phenomenon occurs in SPB stars 
( [Dziembowskj et al.|199"3j|Dupret et al.|2008 1 



Interestingly, the blue edge of the I < 3 instability strip of 
our sdB models with OP opacities corresponds to the observed 
blue edge of the V1093 Her variables at log T e g = 4.47. Thus, 
we have identified that atomic diffusion represents another step 
towards solving the 'blue edge problem', in addition to the use 
of OP opacities and Fe/Ni enhancement suggested by |Jeffery &] 
Saio (2006a b). The closer match between the theoretical and 
observational instability strips is caused by the settling of He 
and the outward diffusion of H. As a result, the envelope's mass 
density decreases and the opacity bump at log T = 5.3 moves 
inwards within the star because of the relation |dT/dr| oc Kp. 
This is towards regions that have a longer thermal relaxation 
time, and as explained before, modes with longer pulsation peri- 
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The results for the first three assumptions are also shown in Fig.|6] 



ods can now be efficiently driven. These are high-order g-modes 
for which no significant radiative damping occurs. For a grid 
that was evolved without H/He-diffusion and with OP opacities 
(lower panels of Fig. |6j, we indeed find the same blue edge as 
|Jeffery & Sato] ( |200 6b), i.e., logr eff = 4.45 for I = 3 modes. 
However, Ni enhancement remains an essential ingredient, as 
we checked for models with only Fe enhancement (see Table 
[TJ. Observations confirm the effect of settling, because most sdB 

stars are He-deficient with log[/V( He)/7V(H)] 2 ( |Saffer et alT] 

[T9941 |Heber & Edelmann|[2t>04l i. However, the effect of He- 
settling is probably less efficient in reality than in our models 
because of counteracting processes such as weak stellar winds 
( |Unglaub & Bues|2001| >, turbulent mixing, and radiative levita- 
tion (although the effect of radiative levitation on He is found to 
be small, see |Michaud & Richer|2008l l. 



4. Conclusions 

The precise details of chemical transitions in subdwarf B stars 
have a considerable influence on the period spacings of high or- 
der g-modes. This is a signature of mode trapping that occurs 
due to discontinuities in the BV-frequency. We have shown that 
the structure of mode trapping depends sensitively on the exact 
shape of the H-profile in the core-envelope transition zone. Thus, 
when constructing ZAEHB models without following the previ- 
ous evolution, the chemical profiles must be modelled carefully. 
Since the action of mode trapping is so sensitive to the shape 
of the chemical transition, it is necessary to include atomic dif- 
fusion during the evolution, i.e., the processes of gravitational 
settling, temperature and concentration diffusion. 

Furthermore, if the sdB star started He-fusion in a run-away 
flash, the inner part of the core is convectively mixed and C is 
produced up to ~7%. The edge of the He-flash-induced convec- 
tive core is accompanied by a chemical transition. We showed 
that this can cause additional mode trapping features resulting in 
a complicated behaviour of the period spacings. Multiple com- 
position gradients lead to simultaneous mode trapping in differ- 
ent regions. The g-mode spectra become very complicated, and 
it will be difficult to directly infer detailed information about 
the composition gradients from observed g-modes. It is there- 
fore crucial to evaluate the importance of different physical pro- 
cesses, and to account for all of them realistically, as we do here. 

Finally, and perhaps most importantly, we report on a possi- 
ble step forward in solving the 'blue edge problem' of the VI 093 
Her instability strip. Our evolutionary sdB models with (i) OP 
opacities, (ii) the inclusion of H/He diffusion, and (iii) a paramet- 
ric Fe/Ni enhancement, show unstable I < 3 g-modes for effec- 
tive temperatures up to 30, 000 K. This is only achieved when all 
three ingredients are included in the computations. It should be 
noted that it remains unclear which Z-values the observed long- 
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Fig. 6. Pulsation periods of stable (green x) and unstable (red +) modes as a function of the effective temperature. The upper panels 
are for sdB models computed with OPAL opacities and the middle panels are for OP opacities. In the lower panels we show the 
results for OP opacities if H/He-diffusion is not included. In all three cases, the abundances of Fe and Ni are increased around 
log T — 5.3 up to a factor 10 in number during the evolution. 



period pulsation modes have. Thus, although we have demon- 
strated the importance of H/He diffusion, solving the blue edge 
problem effectively will require mode identification, in addition 
to a realistic treatment of radiative levitation and other transport 
mechanisms. 
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